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ABSTRACT 



It is proposed that the spectra of so called 'MeV blazars' can be explained in 
terms of previously developed models of the external comptonization of accretion disk 
radiation provided that the structure of the inner and outer parts of the accretion 
disk is different. The electron acceleration is saturated by the inverse Compton losses 
in the inner geometrically thick disk and the outer geometrically thin disk at differ- 
ent maximum energies which causes the appearance of two spectral components, one 
strongly peaked in the MeV energy range and the second one of a power law type 
extending through the GeV energy range. The spectra, computed in terms of such 
a simple geometrical model, are in a good agreement with observations of the MeV 
blazar PKS0208-512. They are consistent with the transient appearance of a strong 
MeV peak, the power law spectrum in the EGRET energy range, and a possible cut-off 
at high energies. 
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1 INTRODUCTION 

The COMPTEL telescope on the board of Compton GRO 
has detected a few blazars which emitted power is concen- 
trated in the MeV energy range (Blom et al. 1995, Bloe- 
men et al. 1995). Two spectral components can be clearly 
identified in such type of sources (e.g. PKS0208-512, Kan- 
bach 1996, Blom et al. 1996). It has been proposed that 
such spectra can originate in a dense relativistic blob 
moving in a jet. In such model the MeV peak is caused 
by the blueshifted annihilation line (Marcowith, Henri 
& Pelletier 1995, Roland & Hermsen 1995, Bottcher & 
Schlickeiser 1996, Skibo, Dermer & Schlickeiser 1997). How- 
ever observations of 3C 273 by the OSSE telescope show 
that in fact the peak can be sometimes present at ~ 0.3 
MeV, which would rather require the redshifted anni- 
hilation line (McNaron-Brawn et al. 1997). Recently it has 
been suggested that these two emission components may be 
created by two different populations of electrons present in a 
jet, which scatter the external quasi-isotropic soft radiation 
(Sikora & Madejski 1996). 

It has been already noted that the spectra of MeV 
blazars could be explained by the emission of relativis- 
tic electrons which comptonize the thin disk radiation 
(Bednarek, Kirk, Mastichiadis 1996a,b (BKM96a,b)). How- 
ever, recently observed very flat transient spectrum in the 
GeV energy range with a cut-off at a few GeV (Kanbach 



1996) is difficult to explain with the assumption on the thin 
disk geometry. In this paper we show that different stages of 
emission of MeV blazars can be explained in terms of pre- 
viously developed model provided that the structure of the 
inner and outer parts of the accretion disk is different. The 
electron acceleration caused, either by a shock propagating 
along the jet or by electric fields induced in the reconnec- 
tion regions, can be saturated by inverse Compton losses. 
However, if the geometry of the accretion disk, surrounding 
central black hole, changes from the thin one (the outer disk 
part) to the thick one (the inner disk part) and the inner 
parts of the disk are relatively hotter, then the acceleration 
of electrons saturates at different energies during jet prop- 
agation through these radiation fields. It results in an ap- 
pearance of two spectral components, the first one strongly 
peaked in the MeV energy range and the second one, of a 
power low type, extending through the GeV energy range. 



2 GAMMA-RAYS FROM SCATTERING OF 
PHOTONS FROM THICK AND THIN DISK 

We consider the standard external Compton scattering 
model for the 7-ray production in active galaxies, in which 
electrons accelerated in the jet lose energy on inverse Comp- 
ton scattering (ICS) process by scattering the accretion disk 
radiation. In general, the accretion disk may be geometri- 
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cally thin (as considered by e.g. Dermer, Schlickeiser & Mas- 
tichiadis 1992, Dermer & Schlickeiser 1993, Coppi, Kartje & 
Konigl 1993, BKM96a,b), or it may be geometrically thick 
(Melia & Konigl 1989, Bednarek 1993, Bednarek & Kirk 
1995 (BK95)). The electrons in turn may be accelerated 
by a relativistic shock moving in the jet (e.g. Dermer & 
Schlickeiser 1993, Sikora, Begelman & Rees 1994, Bland- 
ford & Levinson 1995), or they move almost rectilinearly 
along the jet axis (Mclia & Konigl 1989, Coppi, Kartje & 
Konigl 1993, BK95, BKM96). In this paper we consider the 
simpler case, assuming that the electrons are accelerated 
rectilinearly in the jet in a relatively small localised regions, 
either in electric fields induced in magnetic reconnection 
or by drifting along the surfaces of highly oblique shocks. 
Such shocks may be formed when the jet plasma meets ob- 
stacles e.g., winds of massive stars moving through the jet 
(Bednarek & Protheroe 1997) or debris of close supernova 
explosion (Bednarek 1997). These acceleration regions are 
located along the axis of the jet in our model. 

In the present paper we apply the model for the disk 
geometry considered previously by Melia & Konigl (1989). 
Therefore, we assume that the outer disk, at distances 
greater than rmid from the central black hole, is geomet- 
rically thin and is characterised by the Shakura & Sun- 
yaev (1973) temperature profile. The inner geometrically 
thick disk (disk funnel), extends between rin and r^id- The 
funnel has an opening angle 6. For the outer disk we assume 
the Shakura &: Sunyaev temperature profile of the form 

Tout(r) = rmid(r/rmid)~°''''^, (1) 

where Tmid is the surface temperature of the disk at the dis- 
tance Tmid from the disk centre. For the temperature profile 
of the inner, thick disk, we assume that it has a similar form 
as the outer disk but with different exponent 

Tfunir) = Tmidir/vmidr"- (2) 

We discuss the cases with 5 > 0.75, which means that the 
inner disk is hotter than expected in Shakura & Sunyaev 
model. It is worth to note that this is consistent with the 
observations of the soft X-ray excess in many active galaxies, 
between them 3 C273 which is also detected in 7-rays (Lichti 
et al. 1995). Hence the disk radiation field is defined by four 
parameters: Tmid, rmid, S, and 9. The distribution of the 
electric field along the jet is defined following BKM96b, 

Eiz) = Eoiz/zo)-", (3) 

where Eo is the electric field strength at the distance zo 
along the jet, and the exponent q is a parameter. 

The Lorentz factor of electrons in the jet, 7eg, is de- 
termined by the balance between the energy gains from the 
rectilinear accelerator and the energy losses on ICS in the 
disk radiation as defined in BK95, and BKM96a,b, 

Pics{le^{z))=ceE(z), (4) 

where Pics are the ICS losses of electrons with Lorentz fac- 
tor 7eq at a distance 2, c is the velocity of light, and e 
is the elementary electric charge. Bednarek, Kirk & Mas- 
tichiadis (BKM96a,b) has shown that the jet can be sepa- 
rated in two regions. In the inner part, called the "radiation 
dominated jet", electron energy losses are strong enough 
that they can reach the equilibrium Lorentz factor deter- 
mined by the balance between energy gains and losses. In 
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Figure 1. The ICS 7-ray spectra produced by electrons which 

scatter the radiation of a thick and thin accretion disk. The pa- 
rameters of the disk are following: T^^^ = 5 X lO** K, » ^jd = 2rin, 
Tin = 10"'^'' cm, 5 = 2, and 9 = 8O". The spectra marked by the 
full curves are obtained for Eo = 0.08 V cm^-', /roc = 0.005, and 
different values of a = 0.67, 1., 1.3?>, 1.67 (counting from the top). 
The dotted curves show the spectra for the same value of hec 
but: Eo = 8x 10-3 V cm"! (bottom dotted curves), Eq = 0.8 V 
cm"^ (upper dotted curve), and a = 1.33. 



the outer part, called the "particle dominated jet", the ICS 
losses are too small and electrons are injected into the jet 

with the Lorentz factors corresponding to the maxirrmm po- 
tential drop in the acceleration regions. These two regions 
are separated at Zmax, measured along the jet, which is given 
by the condition (presented in BKM96a,b) , 

) = elcohE{z max 

), (5) 

where Icoh = l3zmax is the coherence length for the acceler- 
ation of electrons (it can be scaled by the distance from the 
black hole with the scaling factor /3 << 1), E{zmax) is the 
electric field in the acceleration region (given by Eq. (3)), 
and me is the electron rest mass. In the radiation domi- 
nated jet, the energy of electrons is directly transferred to 
the 7-ray photons in ICS process. In the particle dominated 
jet the electrons are likely to be isotropized by the random 
component of the magnetic field after emerging from the 
acceleration regions. These electrons lose energy mainly on 
synchrotron process contributing to the lower energy part of 
the blazar spectrum. 

In general, electrons accelerated in the electric fields 
may scatter disk radiation in the Thomson or the Klcin- 
Nishina case. Here we concentrate on the Thomson regime. 
If the scattering occurs in a narrow disk funnel, the Lorentz 
factor of electrons has to be limited to the values deter- 
mined by the condition ^eqi-iksT /nie(?) < 1, where fes is 
the Boltzmann constant and T is given by Eq 2. If the scat- 
tering occurs above the surface of a thin disk, the condition 
derived in BKM96a (see Eq. (5) in that paper) has to be 
fulfilled. 

The spectra of 7-ray photons produced in the radiation 
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Figure 2. The spectra of 7-rays from ICS, produced by electrons propagating through radiation of the inner geometrically thick disk 
and the outer geometrically thin disk, are compared with the observations of MeV blazax PKS 0208-512 (Blom et al. 1996) during May 
- June 93 (on the left) and July 91 - Jan 92 (on the right). The differences between these spectra are only caused by the change in the 
distribution of the electric field profile along the jet defined by the parameter a, which is equal to a = 1.4 for May - June 93 observations, 
and a = 0.67 for July 91 - Jan 92 observations. For the other parameters of the fit see text. 



dominated jot arc obtained by using the following formula 

dN _ Tin r"-°°' diVe dAr(7e,(^)) 

dE^dQdt c J^^ dndt dE-ydt ^' ^ 

where dNe/dCldt is the injection rate of electrons (per sr~^ 
s~^), dN{jeq{z))/dEjdt is the ICS 7-ray spectrum produced 
bj' an electron, with Lorentz factor "feq(z), located at the 
distance z along the jet. This ICS spectrum is obtained by 
integration of Eq. (2.33) in Blumenthal & Gould (1970) over 
the disk radiation field, zo is taken to be equal to rin, and 
Zmax is determined by Eq. (5). 

In Fig. 1 we show an example of the 7-ray spectra cal- 
culated in terms of this model. In order to obtain these 
spectra we fix the disk parameters and investigate their 
dependence on the jet parameters. Note that the spectra 
show two clear components, corresponding to propagation 
of the jet through the radiation of the inner, geometrically 
thick disk (lower energy part of the spectrum) and propaga- 
tion through the radiation of the outer, geometrically thin 
disk (higher energy part of the spectrum). The parameter a, 
which describes the profile of the electric field in the jet, de- 
termines the shape of the spectrum (see full curves in Fig. 1). 
The location of the maximum in the power spectrum is de- 
termined by the strength of the electric field at the base of 
the jet (see dotted curves in Fig. 1). 



3 COMPARISON WITH THE SPECTRUM OF 
PKS0208-512 

PKS0208-512, at the distance corresponding to the redshift 
z = 1.003, has been observed by the Compton GRO a few 
times (e.g. Blom et al. 1996) showing complex and vari- 
able spectrum. In some observations very strong peak was 
detected in 1 3 MeV energy range by the COMPTEL de- 
tector. The power law spectral component extending above 



~ 30 MeV, with the cut-off at a few GeV, has been also de- 
tected by EGRET detector (Kanbach 1996). The spectrum 
is variable with the spectral index changing from 1.67±0.12 
to 2.24 ±0.36 (Blom et al. 1996). We interpret this complex 
spectral shape in terms of our simple one component jet 
model and Melia & Konigl (1989) geometrical disk model. 
It is assumed that the disk parameters has not been changed 
between these two observations. The disk is defined by the 
inner disk radius rtn = lO^** cm, the radius of the thick disk 
rmid = 2ri„, its temperature Tmid = 5 x lO** K at Vmid, 
the exponent 5 = 2, and the opening angle of the thick disk 
13 = 80°. The fits to the spectra of PKS 0208-512 are shown 
in Fig. 2 during different observations (Blom et al. 1996). 
The spectra are obtained for the following parameters of 
the jet: Eo = 0.08 V cm~^ at the base of the jet zo = rtn, 
and Zrec = 0.005. We assumed that the injection rate of elec- 
trons is uniform along the jet. This might be motivated by 
the fact that the fiux of the electrons in the jet should be 
conserved. 

The change in the spectra, during two observations of 
PKS 0208-512, is caused only by changing the distribution 
of the electric field along the jet, defined by the parameter 
a (see Eq. 3). Its value is equal to 1.4, for the observations 
during May - June 93, and to 0.67, for the observations 
during July 91 - Jan 92. The change in these spectra corre- 
sponds to different equilibrium Lorentz factors of electrons 
during propagation in the electric field of the jet. In Fig. 3 we 
present the values of 7oq as a function of the distance z along 
the jet, for two cases discussed above. Note the difference in 
extraction of the electric field energy by electrons propagat- 
ing through the radiation fields of thick and thin disk in both 
cases. This is the reason of anticorrelation of relative power 
emitted at the high energy (above ~ 100 MeV) and the low 
energy parts of the spectrum. In such model the cut-off in the 
spectrum, observed at around a few GeV (Kanbach 1996), is 
caused by inefficient acceleration of electrons in the electric 
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field of tlie particle dominated jet. This interpretation is dif- 
ferent from another possible explanation of such a feature as 
being due to absorption of high energy 7-rays in the soft rsr 
diation (e.g. Bednarek 1993, Sikora, Begelman & Rees 1994, 
BK95, Blandford & Levinson 1995, Pohl 1996). 



4 CONCLUSION 

We propose a simple geometric explanation as responsible 
for the formation of the two component 7-ray spectrum in 
MeV blazars. It is shown that in terms of the standard ex- 
ternal comptonization model of the disk radiation, the two 
components in the 7-ray spectrum can be explained if the 
disk geometry changes from the geometrically thin one, in 
the outer part, to the geometrically thick one, in the inner 
part. Computed 7-ray spectra can well describe the com- 
plex shape of the most prominent MeV blazar PKS0208-512 
(Fig. 2). The model is able to explain the relative changes 
in the strength of the soft MeV component and hard com- 
ponent above ~ 30 MeV. It is also consistent with the high 
energy cut-off in the 7-ray energy range. Since the lower 
7-ray emission is produced closer to the central engine, it 
is likely that the GeV emission, produced further along the 
jet, is delayed in respect to the MeV emission. This model 
can also simply explain the appearance of a bump in the 
blazar spectra at energies below ~ 0.511 MeV (the case of 
3 C273, McNaron-Brawn et al. 1997), which can not be ex- 
plained by the annihilation in the jet model. The peak 
in the 7-ray power spectrum, predicted by our model, is 
only determined by the balance between energy gains and 
energy losses of electrons propagating inside the thick disk 
(see Fig. 1). Therefore, it may appear in principle, at arbi- 
trary energy. 

It is assumed in this paper that the electrons axe accel- 
erated rectilinearly as discussed originally by Bednarek & 
Kirk (BK95). However, we expect that similar two compo- 
nent spectra should appear if electrons, accelerated by a rel- 
ativistic shock propagating along the jet axis, scatter radiar 
tion of the inner and outer disk. The acceleration of electrons 
can be also saturated by the ICS losses and the maximum 
energy of electrons is determined by the balance between ac- 
celeration efHciency and energy losses. Only the geometry of 
the ICS scattering process by electrons distributed isotropi- 
cally in the blob frame becomes more complicated as shown 
by e.g. Dermer & Schlickeiser (1993). 
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